Bath optimization in the Cellular Dynamical Mean Field Theory 
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In the Cellular Dynamical Mean Field Theory (CDMFT), a strongly correlated system is rep- 
resented by a small cluster of correlated sites, coupled to an adjustable bath of uncorrelated sites 
simulating the cluster's environment; the parameters governing the bath are set by a self-consistency 
condition involving the local Green function and the lattice electron dispersion. Solving the cluster 
problem with an exact diagonalization method is only practical for small bath sizes (8 sites). In 
that case the self-consistency condition cannot be exactly satisfied and is replaced by a minimiza- 
tion procedure. There is some freedom in the definition of the 'merit function' to optimize. We 
use Potthoff 's Self-Energy Functional Approach on the one- and two-dimensional Hubbard models 
to gain insight into the best choice for this merit function. We argue that several merit functions 
should be used and preference given to the one that leads to the smallest gradient of the Potthoff 
self-energy functional. We propose a new merit function weighted with the self-energy that seems 
to fit the Mott transition in two dimensions better than other merit functions. 
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I. INTRODUCTION 

The discovery of high-temperature superconductors in 
the late 1980's and the hypothesis that the mechanism of 
superconductivity in these materials is rooted in strong 
electron-electron interactions has stimulated theoretical 
investigation of lattice models of strongly-correlated elec- 
trons, such as the Hubbard model. One of the early suc- 
cesses of this program was a new understanding of the 
metal-insulator transition in the Hubbard model using 
Dynamical Mean-Field Theory (DMFT) 1 . A central hy- 
pothesis behind DMFT, proven exact in the limit of in- 
finite dimension, 2 is that the momentum dependence of 
the electron self-energy S(k, oj) may be neglected; DMFT 
focuses instead on the frequency dependence, which it 
determines approximately within a self-consistent proce- 
dure. This is equivalent to replacing the original Hub- 
bard model by an effective model in which a single site 
- the "impurity" - is embedded in the lattice through 
hybridization with a bath of uncorrelated orbitals. The 
one-particle Green function G(ui) for the correlated site 
then takes the form 



G _1 (w) =u-T(u) -S(w) 



(1) 



where £(w) is the approximate electron self-energy and 
r(w) is the so-called hybridization function that incorpo- 
rates the effect of the uncorrelated bath on the electron 
propagation. The hybridization function is found by an 
interative procedure that involves (i) the solution of the 
impurity model and (ii) self-consistence between the elec- 
tron Green function at the impurity site and the Green 
function G(k, uj) constructed from the self-energy S(w) 
via Dyson's equation (more details on this below). 

The importance of short-range antiferromagnetic fluc- 
tuations in Hubbard models and the possible existence of 
d-wave pairing has motivated the extension of DMFT to 
procedures where not only a site, but a finite cluster of 
sites, is embedded in the full lattice via a hybridization 
function. This is further motivated by the strong momen- 



tum dependence of the self-energy inferred from photoc- 
mission experiments. 3 The Dynamical Cluster Approxi- 
mation (DCA) 4,5 and the Cellular Dynamical Mean-Field 
Theory (CDMFT) 6 are two generalizations of DMFT to 
finite clusters, that take into account short-range spatial 
correlations. Both have revealed antiferromagnetic order 
and d-wave pairing in the Hubbard model. 7-12 They differ 
in that DCA is formulated in reciprocal space, by parti- 
tioning the Brillouin zone into a finite number of patches, 
whereas the CDMFT is formulated in direct space, by 
tiling the lattice into identical clusters with open bound- 
ary conditions. 

Solving the impurity problem in DMFT (or in its clus- 
ter extensions) may be done at finite temperature via 
Quantum Monte Carlo (QMC) or, as proposed in Ref. 13, 
by the exact diagonalization (ED) of an associated An- 
derson Hamiltonian at zero-temperature. The Monte 
Carlo approach has the advantage of simulating an effec- 
tively infinite bath and of providing temperature infor- 
mation. On the other hand, it suffers from the infamous 
fermion sign problem, which makes convergence slow and 
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FIG. 1. (Color online) Examples of clusters (blue circles) 
with finite baths (squares). Cluster (A) is used to treat the 
one-dimensional Hubbard model and (B) the two-dimensional 
Hubbard model. 



the computational requirements important, even with the 
latest continuous-time algorithms free of discretization 
error. 14 It is also very difficult to extract real-frequency 
dynamical information from it. The pros and cons of the 
ED approach complement those of the QMC approach: 
(1) it is usually carried at zero temperature; (2) it pro- 
vides real-frequency information; (3) it does not suffer 
from the sign problem but (4) it is limited to a small, 
discrete bath system. This last characteristic is its most 
serious drawback, and the one we will deal with in this 
paper. 

In Section II we review the CDMFT algorithm and 
explain how the restriction to a finite bath turns the self- 
consistency into an optimization problem whose solution 
depends on a choice of merit function. In Section III we 
review Potthoff 's self-energy functional approach and as- 
sert that it provides the best possible hybridization func- 
tion for the finite-bath problem. In Section IV we show 
the outcome of CDMFT calculations on the one- and two- 
dimensional Hubbard models, using various merit func- 
tions, and compare them with SFA results. We propose 
to select the merit function that minimizes the gradient 
of the Potthoff functional, in the cases where the SFA 
cannot be practically applied. 
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FIG. 2. (Color online) The CDMFT algorithm with an exact 
diagonalization solver. 



exact diagonalization in order to accelerate convergence 
and cut memory costs. 15,17 

The effect of the bath on the electron Green function 
is encapsulated in the so-called hybridization function 



II. THE CELLULAR DYNAMICAL MEAN 
FIELD THEORY 
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(3) 



In CDMFT, the full lattice Hubbard Hamiltonian H is 
replaced by the following cluster Hamiltonian: 



which enters the cluster electron Green function as 



#' = - I>i4c» + &X 



riifnn 



Y M c L a /^ + H - c -) + Yl s^ a U a ^ 



(2) 



Z,fl,(T 



where ci a annihilates an electron of spin uona physical 
site labelled i, and a fla annihilates an electron of spin a 
on a bath orbital labelled \i. The bath is parametrized 
by the energy of each orbital (e^) and by the bath-cluster 
hybridization matrix 9 ifi (we assume spin-independence 
for simplicity; this would not be the case in a treat- 
ment of antiferromagnetism) . Note that 'bath site' is 
a misnomer, as bath orbitals have no physical position 
assigned to them. The clusters and baths used in this 
work are shown on Fig. 1. Our convention is to add to 
the above Hamiltonian a chemical potential term —fxN, 
where N is the total number of electrons in the cluster 
and the bath. 

Note that a more general parametrization of the sys- 
tems illustrated on Fig. 1 is possible 15,16 ; in particular, 
the most general bath made of 8 orbitals and compati- 
ble with the discrete symmetry of the cluster of Fig. IB 
would require 16 parameters. But our purpose is to il- 
lustrate the relationship between the CDMFT-ED proce- 
dure and Potthoff 's self-energy functional approach, and 
accordingly we choose to limit the size of the parame- 
ter set. In practice, cluster symmetries are used in the 
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(4) 



where we hide site and spin indices behind a matrix nota- 
tion for the cluster one-body terms (t', including chemical 
potential), the hybridization function (l~), the self-energy 
(E) and the cluster Green function (C). 

The basic computational task of DMFT approaches is 
to find the best possible embedding of the cluster into the 
original lattice; that is, to find the best possible value of 
the bath parameters. In CDMFT, this is accomplished 
via a self-consistency condition, using the following iter- 
ative algorithm, summarized in Fig. 2: 

1. Start with a guess value of the bath parameters 
(Qiy., £p), that define the hybridization function (3). 

2. Calculate the cluster Green function G'(w) with the 
Exact diagonalization solver and extract the cluster 
self-energy E(w). 

3. Construct the momentum-dependent Green func- 
tion G(k, w) from that self-energy and project it 
back on the cluster via a Fourier transform: 



G(w) = 



1 



kG^(k)-E(«) 



(5) 



where Jj- stands for an integral over the reduced 
Brillouin zone, along with the appropriate factors 
of 2tt. 



4. Calculate the matrix 

5 - 1 (u;) = G- 1 + Z(w) 



(6) 



5. Choose new bath parameters that make the com- 
bination ui — t' — T(cj) as close as possible to Gq ■ 
Since we have a finite number of bath parameters at 
our disposal, this self-consistency condition cannot 
be fulfilled for all frequencies, but only optimized. 
This is done by minimizing the distance function 



sense that its form does not depend on the one-body part 
t of the Hamiltonian, but only on the interaction part. 
This means that it can be evaluated by considering a dif- 
ferent Hamiltonian, H' , that shares the same interaction 
part with the original Hamiltonian H, but has a different 
one-body part. In particular, the cluster Hamiltonian (2) 
falls in that category. This allows us to extract the value 
of the functional F evaluated at the physical self-energy 
of H' and to calculate 



Cl t [Z] =0' + Trln(-G')-Trln(-G) 



(11) 



Y^ W(w) (w - 1' - l~(w) - G l (u)).. (7) whcre G stands for ( G -i _ E )-i and n > is the grand 



over the set of bath parameters. Changing the bath 
parameters at this step does not require a new so- 
lution of the Hamiltonian H' , but merely a recal- 
culation of the hybridization matrix (3) . 

6. Go back to step (2) with the new bath parameters 
obtained from this minimization, and repeat until 
they are converged. 

In practice, the distance function (7) can take various 
forms, depending on the heuristic weight function W{w). 
So far, ad hoc criteria and intuition have been used to 
choose the weight function. Those that are benchmarked 
in this paper are listed in Eq. (15) below. The sum over 
frequencies in Eq. (7) is carried over a set of equally 
spaced Matsubara frequencies, defined by some fictitious 
inverse temperature /?, typically ranging from 20 to 200 
(in units of t~ r ). 



III. THE SFA APPROACH 

There is a more fundamental way to find the best pos- 
sible value of the bath parameters: Potthoff 's self-energy 
functional approach (SFA). 18 Let us summarize this ap- 
proach here: The self-energy E of any system obeys 
a variational principle in the space of all possible self- 
energies: 



ST 



= 



(8) 



where the functional $7 t ,[/[E] is given by 

Cl t ,u[T} = Fu[Y] - Tr lnt-G^ 1 + E) (9) 

Got is the noninteracting Green function and -Fj/[Z] is 
the Legendre transform of the Luttinger-Ward functional 
$[G]: " 



F u [Z\ = Q u [G]-Tr(TG) 



SG 



(10) 



The value of the functional Cl t .u [E] evaluated at the phys- 
ical self-energy is precisely the grand potential CI of the 
system. Moreover, the functional Fu is universal, in the 



potential associated with the cluster Hamiltonian H' . A 
more explicit expression is 



O t [E] = Cl' 



^E 



lndct 



l-(t fc -t')G'(k,u;) (12) 



where t^ is the one-body matrix of the original system 
(here expressed as a matrix over cluster site indices and 
a function over the reduced Brillouin zone) and t' is the 
one-body matrix of the cluster Hamiltonian H' . T is the 
absolute temperature and the sum is carried over Mat- 
subara frequencies (at zero temperature this translates 
into an integral over the imaginary- frequency axis). 

The bath system is of course assumed to be decoupled 
from the cluster in the original Hamiltonian H. With 
a finite bath, the functional Cl t becomes effectively an 
ordinary function of the bath parameters 9^ and e^. An 
objective answer to the question of what are the optimal 
values of the bath parameters is obtained by solving the 
stationary condition dCl t /dh = 0, where h stands for 
any one of the bath parameters. This condition may 
be further explicited as 



E Tr 



g'^H-G-V) 



dh 



= 0. 



(13) 



where G(w) stands for the Brillouin zone averaged Green 
function defined in Eq. (5). Note that the hybridization 
function l~(u;) enters Gq,., , as can be seen from Eq. (4). 
The distance function (7), by contrast, may be recast as 



W(u>) Tr 



g'-Vj-g-V) 



(14) 



The condition (13) does not entail d — 0, since the lat- 
ter can only be satisfied with an infinite number of bath 
parameters, d being a sum of positive-definite contribu- 
tions. The SFA solution is therefore not self-consistent 
(G'(w) ^k G(w)), but it is, in a variational sense, 19 the 
best possible approximation to the original Hamiltonian 
by a cluster-bath system. 

If we think of the expression between brackets in (14) 
as forming a vector in both frequency and site index 
space, then the CDMFT tries to make that vector as 
small as possible, whereas the SFA tries to make it or- 
thogonal to the vector dT.' /dh, or to make it vanish, the 
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FIG. 3. (Color online) Plots illustrating the CDMFT and SFA solutions from the cluster of Fig. 1A, for the one-dimensional 
Hubbard model, all expressed as a function of the chemical potential fi. Panel (A): Electron density n (the exact Lieb & Wu 
result is also shown). Panel (B): SFA gradient associated with the various CDMFT solutions. Panel (C) and (D): Hybridization 
parameters #1,2 and bath energies £1,2 for the weight functions (15) and the SFA solution. Panel (E) and (F): Bath energies 
£1,2 for a sharp cutoff (15a) with various values of uj c , and the corresponding values of the SFA gradient. Unless indicated 
otherwise, the cutoff was set at uj c — 2 and the fictitious inverse temperature at /3 = 100. 



latter being impossible with a finite bath. Whereas a 
perfect solution (d = 0) of the self-consistency condition 
would automatically satisfy Condition (13), the converse 
is not necessarily true. The SFA is not the only func- 
tional formulation from which DMFT can be derived (see 
Ref. 20 for an extensive discussion). However, it has the 
distinction of being based on an exact evaluation of the 
functional, albeit in a restricted space of parameters. 19 



Detailed benchmarks of the SFA approach against the 
exact solution of half-filled the one-dimensional Hubbard 
model have been performed in Ref. 21, for a large vari- 
ety of clusters (with and without baths). In this work 
we focus instead on comparing the SFA approach with 
CDMFT, at, and away from, half-filling. 



IV. RESULTS 



Fig. 1A illustrates a simple bath-cluster system that 
can be used to approximate the one-dimensional Hub- 
bard model: two cluster sites and four bath sites, with 
two hybridization parameters #1,2 and two bath energies 
£1,2 ■ The presence of four bath sites (as opposed to two) 
is required in order to have particle-hole symmetry at 
half-filling. This bath system was used, for instance, in 
Refs 22 and 23. 

From a computational point of view, solving the SFA 
condition (13) is more difficult to carry out than the self- 
consistent CDMFT algorithm of Fig. 2. More instances 
of the impurity solver must be called, and the optimiza- 
tion of the functional (11) with respect to the bath pa- 
rameters must be very carefully done: it requires great 
numerical precision and is prone to instabilities, because 
of the relatively weak dependence of $7 t on the bath pa- 



rameters. Thus, solving the SFA condition (13) can only 
carried out in practice on small systems with few varia- 
tional parameters, but it sheds light on the proper choice 
of weighting function W to be used in the CDMFT self- 
consistent procedure. 

We have tested the following weight functions W(ui) 
against the SFA results: 

W(iu>) = 1 within u £ [0,w c ] (sharp cutoff) (15a) 

W(iuj) = 1/ui (extra weight to low frequencies) (15b) 



W(iu) = Tr|Z 2 (iw)| within w € [0,< 



(15c) 



All of these functions have a finite support between 10 = 
and some cutoff frequency iuj c , and are evaluated on a 
grid of Matsubara frequencies defined by a 'fictitious' 
temperature 1//3 (recall that the ED solver is used strictly 
at zero-temperature in this work). Using a range of fre- 
quencies along a segment parallel to (and slightly above) 
the real axis has also been tried, but gives very unreliable 
results, presumably because the landscape of the distance 
function (7) is much more complicated, as the zeros (and 
poles) of the Green function are located on the real axis. 
Other weight functions W have also been benchmarked, 
for instance by putting greater emphasis on low frequen- 
cies (W(w) = 1/w 2 ), or less emphasis (W(w) = 1/^/u), 
or proportional to | TrE| instead of Tr |Z| 2 . They bring 
nothing qualitatively different than the choices (15) that 
we benchmark here. 

Note that these various distance functions are all in- 
dependent of the choice of basis used for bath or site 
orbitals: The factor multiplying W in (7) is a matrix 
trace, and so is choice (15c). Choice (15b) is motivated 
by the desire to give much more weight to low frequen- 
cies, and Choice (15c) by the desire to give more weight 
to frequencies with a large self-energy The cutoff w c may 
also be chosen so as to weigh more low frequencies. 

Fig. 3 summarizes the benchmarks we have conducted 
on the one-dimensional Hubbard model, using the clus- 
ter illustrated on Fig. 1A, with U = 4 and nearest- 
neighbor hopping t = 1. Panels (C) and (D) show the 
optimal value of the bath parameter #i i2 and £i j2 as a 
function of chemical potential (within each pair, the two 
parameters are of course interchangeable, and so there 
is no point in labelling them separately). The SFA re- 
sult is non-analytic at a value of the chemical poten- 
tial (n c « 1.4) corresponding to the edge of the gap 
in the one-dimensional Hubbard model. In the range 
/i£ [fi c , U/2], the SFA values of e\ + jU and 9\ are con- 
stant, as they should be since the physical state of the 
system is the same for all values of fi within the Mott 
gap (the same applies to £ 2 and 2 ). This supports our 
view that the SFA provides the best possible values of the 
bath parameters. On the other hand, the bath parame- 
ters obtained from CDMFT with a sharp cutoff (15a) or 
a self-energy weight (15c) are analytic at fj, c , even though 
they follow the general trend of the SFA solution. The 
weight function ~ l/u> (15b) leads to solutions that show 
some non-analycity, but depart more from the SFA solu- 
tion than the other two. 
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FIG. 4. (Color online) Gradient of the SFA functional for 
the various CDMFT weight functions (15). Same system as 
Fig. 3, but this time at half-filling, as a function of U . 
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FIG. 5. (Color online) Ground state energy density of the 
one-dimensional Hubbard model, estimated in various ways 
(see text). The exact result is shown for comparison. 



Fig. 3B shows the norm of the gradient of the SFA 
functional (12) evaluated at the CDMFT solutions found 
using the weight functions (15) (this gradient is precisely 
zero at the SFA solution). This can be used as a mea- 
sure of the departure from the SFA solution, even in cases 
where the SFA solution is not known. Again, the weight 
functions (15a) and (15c) appear to be the most sensi- 
ble, while the one that enhance low frequencies (15b) has 
most of the time the largest gradient. 

Fig. 3E shows the CDMFT values of the bath energies 
£i. 2 , using a sharp cutoff (15a), for different values of the 
cutoff frequency ui c . Panel F shows the SFA gradient for 
the same set of data. This clearly shows that u c should 
be small enough, but that lo c — 1 is too small and gives 
non optimal results in the gapped region. As a rule, 
the value uj c = 2 provides the best results. One could 
also display the same type of analysis as a function of 
fictitious temperature /3 _1 . In that case, one can show 
that the value j3 — 100 (i.e., a fictitious temperature at 
1% of the hopping amplitude t) is a good choice, in terms 



of smoothness and ease of calculation; this is the value 
that was used in all other plots of this paper. On the 
other hand, j3 = 20 is definitely too low. 

Finally, Fig. 3A shows the CDMFT values of the elec- 
tron density n, using the weight functions (15), as well as 
the value obtained from the SFA solution and, this time, 
from Licb and Wu's exact solution of the one-dimensional 
Hubbard model 24 . Again, the weight functions (15a) and 
(15c) are closest on average to the SFA solution. Note 
however that the latter does not coincide with the exact 
solution and that the weight function (15b) yields a solu- 
tion that is sometimes closer to the exact solution. The 
SFA solution would move closer to the exact solution if 
either the number of sites or the bath size were increased. 
But we argue that it provides the best solution for the 
cluster and bath used here, and that it should be the 
standard against which the different CDMFT solutions 
are compared. We view the occasional close proximity of 
a CDMFT solution to the exact solution as accidental. In 
particular, the solutions obtained from the low-frequency 
weight function (15b) can be tuned to yield the correct 
value of the critical chemical potential /i c by adjusting /3, 
but that does not mean a convergence towards the exact 
value as {i~ x — > 0. 

Fig. 4 shows the same type of comparison, this time at 
half-filling, as a function of U. Again, the sharp cutoff 
(15a) and the self-energy weight (15c) stand out as the 
best choices. Note that the gradient goes to zero in the 
U/t —> and t/U — > limits, which is natural given that 
quantum cluster methods such as CDMFT become exact 
in these limits. 

Fig. 5 shows estimates of the ground state energy den- 
sity Eq of the ID Hubbard model as a function of density 
n. The exact Lieb & Wu result is shown for comparison, 
as well as the SFA value obtained from the optimal value 
£1 of the functional (12) by the relation Eq = 17 + fin. 
For the sharp cutoff (15a), we provide estimates of the 
ground state energy density obtained in two ways: (1) 
by calculating the average (K + V) , where the average 
(A') of the kinetic energy is calculated from the lattice 
Green function G(ui, k), and the average (V) of the po- 
tential energy is calculated from the ground state double 
occupancy at a cluster site; (2) by calculating the func- 
tional (12) and adding \xn. We conclude from this graph 
that method (2) provides a better estimate of the ground 
state energy than method (1), even though the CDMFT 
solutions are not exact solutions of the variational equa- 
tions (13). The same conclusion is reached with the other 
weight functions (15). 

Let us turn to the two-dimensional Hubbard model, 
again with nearest-neighbor hopping only. The cluster- 
bath system is illustrated on Fig. IB. We will start with a 
discussion of the half-filled system, in which case only two 
bath parameters are necessary because of particle-hole 
symmetry: a hybridization 8 (dashed lines on Fig. IB), 
and a bath energy ±e (positive on bath sites labelled 5 
through 8, negative on the others). The SFA solution for 
this system was obtained in Ref. 25 as a function of U 
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FIG. 6. (Color online) Top: hybridization parameter 9 for 
the half-filled, two-dimensional Hubbard model, with weight 
functions (15a), compared with the SFA result taken from 
Ref. 25. Middle panel: same for the bath energy e. Bottom 
panel: the SFA gradient for the same solutions. 



and revealed a Mott transition with two critical values of 
the Coulomb repulsion (U c i and U C 2). This is illustrated 
by the full black curve on Fig. 6: there is an insulat- 
ing solution at strong U that overlaps with a metallic 
solution as smaller U, and an 'unstable' solution linking 
the two, as one would expect for a first-order transition 
described, e.g., by Landau's theory of phase transitions. 
The CDMFT solutions that are found for the same prob- 
lem are shown on Fig. 6 for three possible weight func- 
tions. The Mott transition is visible through an upturn 
of the bath parameters, but no hysteresis was observed: 
sweeping U upwards or downwards did not make any dif- 
ference. The bottom panel of Fig. 6 shows the SFA gradi- 
ent calculated from the three CDMFT solutions. Overall, 
the weight function (15c) has the lowest gradient, except 
exactly at the transition; paradoxically, this weight func- 
tion also best describes the transition: the minimum of 
9 and the sharp upturn in e occur right between the U c \ 
and U C 2 found in the SFA, and the changes observed near 
the transistions are the sharpest of the three weight func- 
tions used. At the Mott transition, one would naturally 
expect a hysteresis loop from CDMFT solutions, or a 
local increase in the SFA gradient due to shift of the so- 
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FIG. 7. (Color online) Top: d-wave order parameter as a 
function of density in the two-dimensional Hubbard model, 
obtained from CDMFT, with the weight functions (15). Bot- 
tom: the SFA gradient for the same solutions. 



lution from the vicinity of the metallic SFA solution to 
the vicinity of the insulating SFA solution. The solution 
found using the self-energy weight (15c) does precisely 
that. 

Note that the 1/w weight (15b) seems particularly in- 
adequate for the metallic solution, which runs against 
the intuition that the 1/w weight would better describe 
a state with low-energy states like a metal. This is also 
true of the one-dimensional system described in Fig. 3. 
However, this is naturally understood in the context of 
the variational equations (13): in the low-frequency limit, 
the self-energy of a Fermi liquid vanishes whereas that of 
a Mott insulator or pseudo-gapped system is large. Thus, 
if the appropriate weight is to be based somehow on the 
self-energy, the 1/w weight function should be more ap- 
propriate for a Mott insulator, not a metal. 

Finally, we probe d-wave superconductivity in the two- 
dimensional Hubbard model. This is done like in Refs 
9 and 11, using the cluster- bath system illustrated on 
Fig. IB. In this case we used six bath parameters: a pair 
(ei,8i) of bath energy and hybridization for the 'first' 
bath, made up of the orbitals labelled 5 through 8; a 
similar pair (£2,^2) f° r the 'second' bath, made up of the 
orbitals labelled 9 through 12; two pairing parameters 
(d\,d%), multiplying pairing operators di,2 with d-wave 
symmetry, symbolically represented by the dashed curves 
on Fig. IB. The exact expression of d,\ is 



— 05t a 7i — a7t a 5i — a 6 ^a 8 ^ — a8T a 6i + H.c 



(16) 



and d,2 has the corresponding expression for bath sites 9 
to 12. The introduction of the pairing parameter breaks 
the conservation of particle number in the cluster-bath 
system, and anomalous averages may be nonzero, which 
is taken as the signature of superconductivity in the sys- 
tem, in particular, what we call the d-wave order param- 
eter is the ground-state average of the operator 

b = I (cos(fc x ) - cos(/c y ))(c kT c_ ki + c_ k t c i4 + H.cJ 

(17) 

where J. stands for an integral over the original Brillouin 

zone (along with the appropriate factors of 2n) . The av- 
erage (D) can be calculated from the lattice Green func- 
tion G(w,k) obtained from the CDMFT solution. 

Fig. 7A shows the d-wave order parameter for the 
two-dimensional Hubbard model with nearest-neighbor 
hopping t = 1 and on-site repulsion U = 8, as a function 
of electron density. We only show the hole-doped 
side, since this system is particle-hole symmetric. The 
outcome of four CDMFT weight functions is plotted. 
The three weight functions (15) are used with w c = 2, 
and in addition a higher cutoff (w c = 5) is used with 
weight function (15c). At this point this system is 
beyond the reach of the SFA, as it presents considerable 
numerical challenges. However, this does not prevent 
us from computing the SFA gradient once the CDMFT 
solution is found, and this can be used as a heuristic 
measure of the proximity to the unknown SFA solution. 
All CDMFT solutions found have superconductivity, 
and roughly in the same range, so the choice of weight 
function is largely a quantitative, not qualitative, issue. 
The use of a larger cutoff (w c = 5) can be rejected on 
the basis that it displays the largest SFA gradient of the 
set (this agrees again with the conclusions drawn from 
the one-dimensional system). The 1/w weight function 
(15b) seems more appropriate in the underdoped region: 
it has the lowest SFA gradient there and the largest 
SC order parameter; this is also where the spectral gap 
is largest (not shown). On the other hand, the sharp 
cutoff weight (15a) is more adequate in the overdoped 
region, where the SC gap is smaller. This again confirms 
that the 1/w weight should not be used in a metallic of 
low-gap phase. 



V. CONCLUSION 

In the Cellular Dynamical Mean Field Theory, the self- 
consistency condition G(w) = G(w) cannot be exactly 
satisfied when using an exact diagonalization solver, be- 
cause of the small size of the bath. In other words, the 
distance function (7) cannot be made to vanish, but can 
only be minimized. In that case, an ambiguity arises be- 
cause of the arbitrariness in the choice of the weight func- 
tion W(u>), and this ambiguity translates into a variety of 
solutions with sometimes important quantitative differ- 



ences. We argued that Potthoff's Self-Energy Approach 
provides the best possible solution for the CDMFT bath 
parameters; in particular, it is non- analytic at the critical 
value /i c of the chemical potential that separates the Mott 
and metallic phases and the bath hybridizations and en- 
ergies are independent of /z within the Mott gap; also, the 
SFA yields a first-order Mott transition as a function of U 
in two dimensions, which is not seen in the CDMFT solu- 
tion of the same bath system. We benchmarked different 
weight functions against the solution found in Potthoff's 
Self-Energy Approach, and argued that the best choices 
are provided by weight functions that have the lowest 
SFA gradient. Weight functions that promote low fre- 
quencies in an 'exaggerated' way, e.g. as 1/w (15b), are 
less adequate for phases with no or weak spectral gap, be- 
cause the self-energy is small at low frequencies in those 



phases. We proposed a weight function proportional to 
the self-energy squared (15c), inspired by the role played 
by the self-energy in the SFA variational condition (13). 
This new weight function is the most successful at de- 
scribing the [/-driven Mott transition in two dimensions. 
We also pointed out that the fictitious temperature used 
in evaluating the distance function should be sufficiently 
small, j3 = 100/t being a good rule-of-thumb value. 
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